Magnetic inversion symmetry breaking and ferroelectricity in TbMnOs 



IT) 

o 
o 

c 

m 

(N 



i 

o 
o 



> 

o 

(N 
\D 
O 

in 
o 

03 



c 

o 
o 



X 
S3 



M. Kenzelmann 1 ' 2 ' 3 , A. B. Harris 4 , S. Jonas 3 , C. Broholm 3 ' 5 , J. Schefer 2 , 

S. B. Kim 6 , C. L. Zhang 6 , S.-W. Cheong 6 , 0. P. Vajk 5 and J. W. Lynn 5 
(1) Laboratory for Solid State Physics, ETH Honggerberg, CH-8093 Zurich, Switzerland 
(2) Laboratory for Neutron Scattering, ETH Zurich & Paul Scherrer Institute, CH-5232 Villigen, Switzerland 
(3) Department of Physics and Astronomy, Johns Hopkins University, Baltimore, Maryland 21218 
(4) Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104S 
(5) NIST Center for Neutron Research, National Institute of Standards and Technology, Gaithersburg, Maryland 20899 
(6) Department of Physics & Astronomy, Rutgers University, 136 Frelinghuysen Rd., Piscataway, NJ 08854 

(Dated: February 2, 2008) 

TbMnOa is an orthorhombic insulator where incommensurate spin order for temperature Tn < 
41 K is accompanied by ferroelectric order for T < 28 K. To understand this, we establish the 
magnetic structure above and below the ferroelectric transition using neutron diffraction. In the 
paraelectric phase, the spin structure is incommensurate and longitudinally-modulated. In the 
ferroelectric phase, however, there is a transverse incommensurate spiral. We show that the spiral 
breaks spatial inversion symmetry and can account for magnetoelectricity in TbMnOa. 

PACS numbers: 75.25.+Z, 75.10.Jm, 75.40.Gb 



The coexistence of antiferromagnetism and ferroelec- 
tricity in solid materials is rare, and much rarer still is 
a strong coupling between these two order parameters 
0,0. In non-magnetic perovskites like BaTiC>3, ferro- 
electricity is driven by a hybridization of empty d or- 
bitals with occupied p orbitals of the octahedrally coor- 
dinated oxygen ions |3( ■ This mechanism requires empty 
d orbitals and thus cannot lead to magnetic ferroelec- 
tric materials. In materials such as BiMnC>3, lone s 2 
electron pairs can lower their energy by hybridizing with 
empty p orbitals |4J. While it leads to coexistence of 
magnetic order with electric polarization at low tempera- 
tures, the very different ordering temperatures show that 
the two order parameter lower the symmetry of the sys- 
tems through distinctly different collective effects. 

TbMnC>3 is an antiferromagnet which contains Mn 3+ 
ions with occupied d orbitals and no lone s 2 cation 
So^as for a number of recently discovered multiferroics 
(a, [D, la] , neither of the mechanisms described above can 
explain the coexistence of magnetic and electric order. In 
these materials, a magnetic field of a few Tesla can switch 
the direction of the electric polarization |5j - proof of a 
strong direct coupling between the magnetic and elec- 
tric polarization. Rare earth manganese oxides show a 
plethora of exciting phenomena which arise from com- 
peting interactions. LaMnC>3 is the parent compound 
to a series of materials featuring colossal magnetoresis- 
tance, and shows orbital ordering of its e g orbital, giving 
rise to layered antiferromagnetic ordering. With decreas- 
ing rare-earth size, there is a tendency towards incom- 
mensurate magnetic order. Ferroelectric ground states 
were recently predicted for doped manganite of the type 
Ri_ x Ca x MnC>3, but they have not yet been observed 
The coexistence and strong coupling of ferroelectricity 
and antiferromagnetism in TbMnC>3 suggests the pres- 
ence of a non-conventional coupling mechanism involv- 



ing competing spin interactions, and it is thus of both 
practical and fundamental importance. 

To develop a microscopic theory of the new coupling 
mechanism, unambiguous determination of the symme- 
try of the magnetic order is essential. In this Letter, we 
present neutron diffraction measurements of orthorhom- 
bic TbMnC>3, which determine the magnetic ground 
states and the phase diagram as a function of temper- 
ature and a magnetic field H||a. We find that the ferro- 
electric phase transition coincides with a magnetic tran- 
sition from a longitudinal incommensurate structure to 
an incommensurate spiral structure that breaks spatial 
inversion symmetry. We show that a recent theory pro- 
posed for axial non-axial parity breaking [j"ollll| predicts 
the observed orientation of the ferroelectric polarization 
based on the symmetry of the magnetic structure that 
we report. 

TbMnOs crystals were grown using an optical float- 
ing zone furnace. A measurement of the temperature 
dependence of the dielectric constant of the two sam- 
ples which we studied confirmed ferroelectricity below 
T=26 K and T=28 K for sample 1 and 2, respectively. 
The discrepancy may result from slightly differing oxy- 
gen partial pressure during annealing. The space group 
is #62, and in the Pbnm setting the lattice parameters 
are a=5.3A, b=5. 86A and c=7A9A. The measurements 
were performed on two single-crystals weighing 40 mg 
and 220 mg, respectively, with the BT2 and SPINS spec- 
trometers at NIST, and the TriCS 4-circle diffractometer 
at PSI. 

The wave-vector dependence of the diffraction in- 
tensity along the (0, k, 1) direction of reciprocal space 
(Figs.n an( i|U illustrates the temperature dependence of 
magnetic order in TbMnOs. Mn 3+ spins develop long- 
range order at Tjv=41 K Immediately below the 
transition temperature, magnetic Bragg peaks are ob- 
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FIG. 1: Diffraction intensity aiong the (0, fc, 1) direction for 
different temperatures (with a factor 10 offset between the 
sub-Tjv data sets), showing fundamental and higher order 
Fourier components of the magnetic ordering. The 3t and 9t- 
order peaks associated with the Tb order occur at k = 0.725 
and 0.175, respectively, and are difficult to identify unambigu- 
ously due to their proximity to strong peaks. 

served at (0, 1 — q, 1) and (0, q, 1) positions. At 35 K the 
magnetic order is described by a single Fourier compo- 
nent associated with a wave- vector (0, q, 0) with q = 0.27, 
and no apparent higher order Fourier components. Upon 
cooling below T=28 K, third-order magnetic Bragg peaks 
(0, 1 — 3q, 0) appear (Fig.^l, indicating step-like modula- 
tion of the magnetic moment. The inset to Fig. [2b shows 
that the wave vector continues to vary with temperature 
throughout the ferroelectric phase. 

Below T—7 K, additional magnetic peaks appear that 
are associated with Tb moments and indicate that their 
interactions favor a different ordering wave-vector than 
for the Mn moments. At T—2 K, there is a strong peak 
at (0, t, 1) with t — 0.425, and several higher-order peaks 
as indicated in Fig. ^ The many strong odd high-order 
reflections indicate that the order associated with the 
(0, t, 0) wave-vector is strongly distorted (bunched struc- 
ture) as expected for anisotropic rare-earth magnets. The 
correlation length and the incommensuration of the Tb 
order at low-temperatures is sample dependent, as is of- 
ten found in systems with phase transitions between in- 
commensurate structures. The correlation length along 
the b-axis was 58(20) and 280(20)A, and t was 0.41 and 
0.425 in sample 1 and 2, respectively. 

The magnetic structures at T=15 K, and 35 K were 
determined from up to 922 first-order magnetic Bragg 
peaks each. Only first order peaks were included in 
the refinement as these are sufficient to determine the 
symmetry of the magnetic structure. Representational 
analysis was used to find the irreducible representations 
that describe the magnetic structures. The structure 
at T=35 K in the high-temperature incommensurate 
(HTI) phase can be described by a single irreducible 
representation. The best fit with x 2 — 0-86 was ob- 
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FIG. 2: (a) T-dependence of the integrated intensity of the 
(0, q, 1) Bragg peak and the associated third order peak at 
(0, 1 — 3q, 1). (b) T-dependence of the wave- vector q, corrected 
for the T-dependence of the b lattice parameter, showing the 
absence of a lock-in transition at T=28 K. 

tained for a structure described by representation T3, 
and the excellent agreement between the model and 
the data is shown in Fig. . The magnetic structure 

(0.0(8), 2.90(5), 0.0(5)) fJ, B and 



is described by m^ 1 " 



-.Tb 



(0, 0, 0.0(4)) fiB where the subscript denotes the 
irreducible representation indicated in Fig.^t. The mag- 
netic structure is longitudinally-modulated with moment 
along the b-axis, as illustrated in Fig. 0Ji an d consistent 
with an earlier study 01 ■ The absence of observable 
higher order peaks indicates that the magnetic structure 
at T=35 K is sinusoidally modulated. 

Two irreducible representations are required to de- 
scribe the magnetic structure at T=15 K in the low- 
temperature incommensurate (LTI) phase. We found 
best agreement with x 2 — 2.19 for magnetic ordering 
involving T2 and T3, as shown in Fig. [3Jd. Fits us- 
ing the Ti and T3, or the T3 and T4 representation 
pairs led to x 2 — 14-5 and higher, and can thus be ex- 
cluded. Neglecting higher order reflections, the magnetic 
structure is given by mf n = (0.0(5), 3.9(1), 0.0(7)) /i B , 
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m 

and mj b 



(0,0, 0(1)) mf n = (0.0(1), 0.0(8), 2.8(1)) Ms 
(1.2(1), 0(1), 0)/xg. The experiment was not 
sensitive to the phase between the y and ^-component of 
the Mn moment. From the size of the moment, however, 
we deduce that the Mn moments form an elliptical spiral. 
The data did not favor a phase difference between the 
Tb and Mn moments, so these phases remain undeter- 
mined. Symmetry splits the Tb moments into two orbits 
which representation theory normally treats as indepen- 
dent. However, as suggested by Landau theory 0, we 
took these two Tb amplitudes to be identical. The phase 
between the two orbits was found to be to 1.3(3)7r. The 
greatly improved fit is evidence that the Tb sublattice 
carries significant magnetization in the LTI phase, pre- 
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FIG. 3: Integrated intensities of magnetic Bragg peaks mea- 
sured at T=15 K and T=35 K using TriCS are compared to 
various magnetic structures. Peaks are sorted by decreasing 
measured intensity, (a) The T—35 K structure is described by 
a single irreducible representation, T3, and the moments point 
along the b-axis. (b) The T—15 K magnetic structure is a spi- 
ral described by a i/-component of F3 and an x-component of 
T2 with a 7r/2 phase shift. 

sumably as a consequence of the exchange field from the 
ordered Mn sublattice. 

Fig. [SJa-b) shows the field dependence of the (0, q, 1) 
magnetic Bragg reflection, which arises from the Mn spin 
spiral. Both the position and the intensity are field inde- 
pendent to within errorbar - evidence that the structure 
remains a spiral up to at least 6 T. Our calculations show 
that the intensity should drop by 7% if the z-component 
of T2 were extinguished. In contrast, no decrease is ob- 
served to within an error bar of 2% between and 6 T. 

The (0, 1 — q, 1) Bragg reflection shown in Fig.|SjL arises 
from Mn T3 magnetization along the a-axis and from 
Tb T2 magnetization along the a-axis. Because the x- 
component of the Mn moment is small, the (0, 1 — q, 1) 
Bragg reflection is particularly sensitive to Tb order. For 
T < 28 K the (0, 1 — q, 1) intensity is suppressed by a 
field H||a confirming that the modulated Tb moment is 
oriented along that direction. Below the Tb ordering 
temperature, the field dependent magnetic Bragg inten- 
sity has a finite-field maximum (Fig. indicating a 
spin-flop transition. 

We collected 51 magnetic Bragg peaks at T=4 K and 
H=4 T along a to determine the magnetic structure 
at low temperatures above the critical field for (0, t, 0) 
Tb order (Fig. [5J;) . The magnetic structure can be de- 
scribed by T2 and T3 with x 2 = 3.81 or by Ti and T3 
with x 2 = 4.19. Since a field along the a-direction dis- 
favors antiparallel spin alignment in the same direction 
as in the Ti-Fs structure, we infer that the structure 
is given by mf n = (0.3(4), 4.7(3), 0.0(5)) Ms, mj b = 



FIG. 4: Schematic of the magnetic structure at (a) T=35 K 
and (b) T—15 K, projected onto the b-c plane. Filled ar- 
rows indicate direction and magnitude of Mn moments. The 
longitudinally-modulated phase (a) respects inversion sym- 
metry along the c axis, but the spiral phase (b) violates it, 
allowing an electric polarization (unfilled arrow), (c) Irre- 
ducible representation of the group Gk for the incommensu- 
rate magnetic structure with k = (0, q, 0). (d) Phase diagram 
as a function of temperature and field applied along the a-axis. 
Full circles indicate second order phase transitions. Open cir- 
cles indicate the characteristic field for reduction of magnetic 
(0, 1 — q, 1) Bragg scattering from Tb moments by 50% from 
its zero-field intensity. 

(0,0, 0.0(3)) Ms, mf n = (0.0(2), 0.0(4), 3.0(3)) fx B and 
mj b = (0.3(2), 0.0(4), 0)^ B . This result suggests that 
the spin spiral structure is more stable for fields along 
the a-axis than for fields along the b-axis 

Harris et al. 0, recently showed that insula- 
tors with axial-non-axial parity breaking magnetic phase 
transitions must also be electrically polarized. Given the 
magnetic structure and the temperature dependence of 
the magnetic order parameter, the theory predicts the di- 
rection and temperature dependence of the electric polar- 
ization resulting from a symmetry allowed trilinear cou- 
pling term. In the following we show that this theory 
correctly accounts for the direction of the electric polar- 
ization in the LTI phase of TbMnOa, and the absence 
of electric polarization in the HTI phase. The trilinear 
magnetoelectric coupling term in the Landau free energy 
expansion is written as V = ^2 UVJ a UV7 a u (k)a v (— k)P 1 . 
Here u u (k) is the magnetic order parameter of irre- 
ducible representation r„, P 7 is the electric polariza- 
tion along the the 7 crystallographic direction and a UV7 
parametrizes the strength of the interaction between the 
electric and magnetic order parameters. In the HTI 
phase, the magnetic order is described by only one irre- 
ducible representation, T3 and due to the high symmetry 
of the Mn moments, it is possible to define a u such that 
under inversion a u (k) — > a u (k)*. The trilinear coupling 
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FIG. 5: Field dependence of magnetic Bragg scattering from 
TbMnOs. (a) and (b) show data for the incommensurate 
peaks that occur for T < 41 K. (a) The (0, q, 1) peak that is 
mostly sensitive to staggered magnetization on Mn sites and 
the (0, 1 — q, 1) peak that is sensitive to staggered magnetiza- 
tion on Tb sites, (c) and (d) show data for the incommensu- 
rate peaks that develop for T < 7 K. 

thus consists only of terms such as V = J2 7 ctyW3{k)\ 2 Py 
Since interactions in a Landau expansion must have the 
symmetry of the paramagnetic phase, this interaction 
must be invariant under inversion. This requires that 
a 7 vanishes, so that P 7 = in the HTI phase. The con- 
clusion remains valid for Tb order with no restriction on 
the phase between its two orbits. 

For the LTI phase, which is described by r 2 and 
T%, there are, however, additional terms such as V = 
'^2 1 b 1 a2{k)<j^(~k)P 1 + c.c, where c. c. denotes the 
complex conjugate. For V to be an invariant, P 7 must 
transform as the product of r 2 and T3. That is, the elec- 
tric polarization must be even under 1 and m yz , and odd 
under 2 y and m xy . This condition can only be satisfied 
for an electric polarization along the c-axis. Previous 
dielectric experiments have shown that the electric po- 
larization that develops below the transition to the LTI 
phase is indeed oriented along the c-axis. 

The coupling of the magnetic order parameter to elec- 
tric polarization in TbMnC>3 Q is similar to that in 
M3V2O8, which adopts two different incommensurate 
magnetic structures |l4| . one of them ferroelectric and 
described by two irreducible representations. The correct 
prediction of the electric polarization for both TbMnC>3 
and Ni 3 V 2 8 (Ref.O) suggests that magnetoelectricity 
resulting from a trilinear magnetoelectric coupling term 
may be commonplace in insulating transition metal ox- 
ides with non-collinear incommensurate structures. Ac- 
cordingly we find less compelling the suggestion [5] that 
the appearance of ferroelectricity is associated with an 
incommensurate to commensurate phase transition. In- 
deed, there is no evidence of a lock-in transition (see 
Fig. |2Jd) and it therefor appears to be irrelevant from the 



point of view of ferroelectricity whether the modulated 
magnetic order is commensurate or truly incommensu- 
rate. 

In summary, we have determined the magnetic struc- 
ture of the paraelectric and ferroelectric phases of 
TbMnC>3. We showed that the paraelectric, magneti- 
cally incommensurate phase has sinusoidally-modulatcd 
collincar magnetic order that does not break inversion 
symmetry The ferroelectric phase, however, has non- 
collinear incommensurate magnetic order described by 
two irreducible representations, which explicitly breaks 
inversion symmetry and thus gives rise to electric polar- 
ization. The qualitative aspects of magnetoelectric ef- 
fects in TbMnC>3 appear to be accounted for by a tri- 
linear coupling term in a Landau free energy expansion 
as proposed by Harris et al. 0, [uj. Understanding 
the magnitude of the effect will require experimental as 
well as theoretical work to track lattice, charge, and or- 
bital degrees of freedom through axial-non-axial parity 
breaking phase transition in insulators such as TbMnOs. 
Apart from the fundamental challenge, improved under- 
standing of magnetoelectricity in these systems may help 
to produce materials for room-temperature applications. 

We thank A. Aharony, O. Entin-Wohlman, and A. 
Ramirez for helpful discussions. This work was supported 
by the Swiss National Science Foundation under Con- 
tract No. PP002-102831. Work at Johns Hopkins Uni- 
versity was supported by the DoE through DE-FG02- 
02ER45983. Work at University of Pennsylvania was 
supported by the U.S. -Israel Binational Science Foun- 
dation under Grant number 2000073. Work at Rutgers 
University was supported by the NSF-DMR-MRSEC-00- 
080008. This work is based on experiments performed at 
the Swiss spallation neutron source SINQ, Paul Scherrer 
Institute, Villigen, Switzerland. The work at SPINS is 
based upon activities supported by the National Science 
Foundation under Agreement No. DMR-9986442. 



[1] 
[2] 
[3; 

K 

[5] 

[6] 
[7] 
[8 
[9 
[10 

[11 
[12 
[13 
[14 



A. Smolenskii et al, Sov. Phys. Usp. 25, 475 (1982). 
Schmid, Ferroelectrics 62, 317 (1994). 
E. Cohen, Nature 358, 136 (1992). 
. Atanasov et al, J. Phys. Chem. A 105, 5450 (2001). 
Kimura et al, Nature 426, 55 (2003). 
Lottermoser et al, Nature 430, 541 (2004). 
Hur et al, Nature 429, 392 (2004). 
Kobayashi et al, J. Phys. Soc. Jpn 73, 1593 (2004). 
V. Efremov et al, Nature Materials 3, 853 (2004). 
B. Harris et al, to be published. 
Lawes et al, cond-mat/0503385 
Quezel et al, Physica B 86-88, 916 (1977). 
Kajimoto et at, Phys. Rev. B 70, 012401 (2004). 
Lawes et al, Phys. Rev. Lett. 93, 247201 (2004). 



